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Abstract 



We consider the multidimensional Euler-Poisson equations with non-zero heat 
conduction, which consist of a coupled hyperbolic-parabolic-elliptic system of bal- 
ance laws. We make a deep analysis on the coupling effects and establish a local 
well-posedness of classical solutions to the Cauchy problem pertaining to data in 
the critical Besov space. Proof mainly relies on a standard iteration argument. To 
achieve it, a new Moser-type inequality is developed by the Bony' decomposition. 

Keywords: local well-posedness, Euler-Poisson equations, classical solutions, Besov spaces 
AMS subject classification: 35M10; 35Q35; 76X05 

1 Introduction and main results 

The ongoing miniaturization of semiconductor devices, some high field phenomena such as hot 
electron effects, impact ionization and heat generation appear inside the devices. The traditional 
drift-diffusion model employed for numerical simulation does not provide an adequate description 
of these effects. Consequently, the hydro dynamical model for semiconductors was introduced, 
which can be derived from the Boltzmann equation by moment method based on the shifted 
Maxwellian ansatz for the equilibrium phase space distribution. Precisely, the hydrodynamical 
model takes the form of the following compressible Euler-Poisson equations (see, e.g., p6j): 




w-w 




* E-mail: jiangxu_79@yahoo.coin.cn, jiangxu_79@nuaa.edu.cn 



for {t, x) e [0, +00) xR^{N = 2, 3). Here, n{t, x) > denotes the electron density, u(t, x) e 
electron velocity and W{t,x) energy density. = x) represents the electrostatic potential 
generated by the Coulomb force from electrons and background ions. P = riT is the pressure 
of electron fluid where 7~(t, x) is the temperature of electrons. The energy density W satisfies 
W = — h T^il > 1) and W = is the ambient device energy, where 71 > is a given 
ambient device temperature, n > is the doping profile which stands for the density of fixed, 
positively charged background ions. The scaled coefficient Tp, Tyj and A are the momentum 
relaxation-time, energy relaxation-time and the Debye length, respectively. The coefficient k is 
the heat conductivity, which generally depends on the electron density and temperature. For 
the sake of simplicity, we assume it to be one constant. 

The full hydrodynamical model (jl.ip for the balance laws of the density, velocity, temperature 
and the electric potential consists of a quasi-linear hyperbolic-parabolic-elliptic system, which 
contains the damping relaxation, heat conduction and electric dissipation. The interaction of 
these special effects makes it complicated to understand the qualitative behavior of solutions, 
however, many efforts were made by various authors, see [U [H [3l [71 [HI [HI [151 [Ml [13 [13 [HI 1221 
[23l [25l [27l [28l [29l [33l [35] and the references therein, for issues of well-posedness of steady-state 
solutions or classical, large time behavior and singular limit problems. 

In this paper, we are concerned with the well-posedness of classical solutions starting with 
smooth initial data under the coupled effects. In the one space dimension, Chen, Jerome, 
and Zhang [7J first considered the initial boundary problem of (jl.ip and established the local 
existence of smooth solutions. Furthermore, they showed that smoothness in local solutions can 
be extended globally in time for smooth initial data near a constant state. This result indicated 
that the relaxation effect could prevent the development of shock waves for the case of smooth 
initial data with small oscillation. Hsiao and Wang |19j considered the corresponding non- 
constant steady state solutions to (11. ip and it was shown that the solutions were exponentially 
locally asymptotically stable. For the Cauchy problem of (jl.ip with large smooth initial data, it 
was proved in [28j the solution generally develops a singularity, shock waves, and hence no global 
classical solution exists, which is due to the strong hyperbolicity, even the damping relaxation 
and the heat conduction (parabolicity) can't prevent the formation of singularity. The effect of 
the Poisson coupling (ellipticity) is smoothing and it decisively affects the stationary states of 
the Euler- Poisson equations (jl.ip , see [HI [H] ■ 

Physically, it is more important and more interesting to study (jl.ip in several space di- 
mensions where were expected to get some similar results as the one-dimension case. Hsiao, 
Jiang and Zhang |17] first studied the Cauchy- Neumann problem of (jl.ip . Using the classical 
energy approach, they established the global exponential stability of small smooth solutions 
near the constant equilibrium. Subsequently, Li [22j extended their results to the non-constant 
equilibrium. 

Recently, we started a program to investigate the hydrodynamic model for semiconductor 
from the point of view of Fourier analysis, which is more careful and refined manner. For 
instance, the method enables us to understand the Poisson coupling effect well, which plays a 
key role in the low frequency of density. Such a fact explains the global exponential stability of 
small smooth solutions in essential. Up to now, we have achieved some results in this direction, 
see [121 [301 [3TI [32I I34j . In these works, we focused on the case of k = mainly and the Euler- 
Poisson equations (jl.ip can be reduced to the pure hyperbolic form of the balance law with a 
non-local source term by virtue of the Green's formulation. However, the case of k 7^ is not 
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the trivial one based on the following considerations: 

(1) The Euler-Poisson equations (jl.ip has not scaling invariance, pertinent to the compressible 
Navier-Stokes equations in [SJ [10]. Thus, we are going to choose the non-homogeneous 
Besov spaces Bp^l^il < p < oo) as the basic functional setting (in x), which are the critical 
spaces embedding in the space of Lipschitz functions; 

(2) For the critical case of regularity index, the classical existence theory for generally hy- 
perbolic systems established by Kato and Majda [21: i24j fails. Iftimie [20] first gave the 
contribution for general hyperbolic systems, however, due to the complicated coupling, 
the result of Iftimie can not be applied directly; 

(3) The inequality in Proposition 12.61 giving parabolic regularity in the framework of Besov 
spaces depends on given time T (except for the case of a = 1 and ai = oo ), which may 
preclude from proving global existence results (even if small data). 

Therefore, as the first step, we establish a local existence result of Cauchy problem pertaining 
to data in the critical Besov spaces for the Euler-Poisson equations (jl.ip with ac 7^ 0. For this 
purpose, the initial conditions for n, u and T, and a boundary condition for are equipped: 

(n,u,r)(x,0) = (no,uo,ro)(x), x G M^, (1.2) 
lim ^(t,x) =0, a. e. t> 0, (1.3) 

\x\—^ + OD 

where the homogeneous boundary condition for means that the semiconductor device is in 
equilibrium at infinity. 

The local existence of a solution stems from the standard iterative method. In comparison 
with that in |10) . there are some differences in the proof of local existence. First, the estimate of 
density n does not follow from the estimates for the transport equation in Besov spaces directly, 
since the velocity u in the momentum equations has not higher regularity. Actually, to obtain 
the desired frequency-localization estimate of (n,u), we introduce a function change to reduce 
(|l.ip to a part symmetric hyperbolic form, and take full advantage of linear hyperbolic theory in 
the framework of Besov spaces and hyperbolic energy method for dyadic blocks. Consequently, 
this restricts us to the space case of p = 2 in (1). Second, we show that the approximate 
solution sequence is a Cauchy sequence in some norm to prove the convergence rather than using 
compactness arguments. In the meantime, in order to overcome the difficulty arising from the 
heat conduction term, we develop a more general version of the classical Moser-type inequality in 
Proposition l2.3l the reader is refer to the Appendix. According to the new Moser-type inequality, 
the heat conduction term can be estimated ultimately, for details, see ()3.34p - ()3.35p . Finally, by 
a careful analysis on the coupling effects, the uniqueness of classical solutions is shown in the 
appropriately larger spaces. 

Through this paper, the regularity index a = 1 + N/2. Our main result is stated as follows. 

Theorem 1.1. Let n^Ti > be the constant reference density and temperature. Suppose that 
no - n, uo, V<?(-, 0) € ^^^(M^) with uq > and To - Tl & B^^'^iR^). Then there exists a time 
Ti > such that 
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(i) Existence: the system U.l\) - nT3\) has a solution (n, u, T, V<P) belongs to 

(n,u,r, V^) G C^([0,Ti] X M^) with n > for all t £ [0,Ti]. 
Furthermore, the solution (n,u, T, V^P) satisfies 

(n-n,u,V^) GCT,(i32",i(K^)) x (CT,(i?2^(K^))) x (CT,(i?2^(K^))) 

and 

T-Tl€Ct,{B^^H^''))- 

(ii) Uniqueness: the solution (n,u, 7", V<^) is unique in the spaces 

L¥AB^') X (l^,{B^')Y X L},^{Bl,) X (L^^{B^')y, 

In addition, there exists a constant Cq > depending only on n,TL, N,j such that 

|I(n-n,u,V#)||£^^(^^._^) + ||r-rL||£oo(5.+i) < CoM, (1.4) 

where M := ||(no - n, uq, V#(-, 0)) + \\To-Tl\\b-+^ and V<P{-,0) := VA'i(no -n). 
Remark 1.1. The symbol VA~^ means 

VA-V= / V,,G{x-y)f{y)dy, 
where G{x,y) is a solution to AxG{x,y) = 6{x — y) with x,y £ R^. 

Remark 1.2. Unlike [10], no smallness condition on the initial density is required, hence the 
local well-posedness in critical spaces holds for any initial density bounded away from zero. Let 
us mention that the heat conductivity k in the proof is assumed to be large appropriately, see 
the following inequalities (13.2ip and (13.36p . 

Remark 1.3. The local existence results in the framework of Sobolev spaces with higher 
regularity (s > 1 + N/2) were obtained by the standard contraction mapping principle, see, 
e.g., [7|. Theorem 11.11 deals with the limit case of regularity (a = 1 + N/2), which is a natural 
generalization of their results. To the best of our knowledge, the global well-posedness and large- 
time behavior of (II. ip in critical spaces still remain unsolved, since the inequality of parabolic 
regularity in Proposition 12.61 depends on the time T, which are under current consideration. 

The rest of this paper unfolds as follows. In Section 2, we briefly review the Littlewood-Paley 
decomposition theory and the characterization of Besov spaces and Chemin-Lerner's spaces. 
Section 3 is dedicated to the proof of the local well-posedness of classical solutions in critical 
spaces. Finally, the paper ends with an appendix, where we prove a Moser-type inequality with 
aid of the Bony's composition. 

Notations. Throughout this paper, C > is a harmless constant. Denote by C([0,T],X) 
(resp., C^{[0,T], X)) the space of continuous (resp., continuously differentiable) functions on 
[0, T] with values in a Banach space X. For simplicity, the notation || (a, b, c, d) \\x means ||a||x + 
||b||x + ||c||x + |jc^||x, where a, b,c,d € X. We shall omit the space dependence, since all functional 
spaces are considered in M^. Moreover, the integral f^js, fdx is labeled as f f without any 
ambiguity. 
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2 Tools 



The proofs of most of the results presented in this paper require a dyadic decomposition of 
Fourier variable. Let us recall briefly the Littlewood-Paley decomposition theory and the char- 
acterization of Besov spaces and Chemin-Lerner's spaces, see for instance |5l[8j for more details. 

Let ((/?, x) be a couple of smooth functions valued in [0, 1] such that ip is supported in the 
shell C(0, f , I) = U e K^ll < lel < §}, X is supported in the ball B(0, |) = G < |} 

and 

oo 

Let S' be the dual space of the Schwartz class S. For / G S', the nonhomogeneous dyadic blocks 
are defined as follows: 

A_i/ := x{D)f = oj*f with oj = 

A J := ip{2-'^D)f = 2"'^ I u{2'^y)f{x - y)dy with u = J-^ V, if > 0, 

where * the convolution operator and J-"^^ the inverse Fourier transform. The nonhomogeneous 
Littlewood-Paley decomposition is 

f=Y,\f V/ G S'. 

q>-l 



Define the low frequency cut-off by 



p<q-l 



Of course, Sq/ = A_i/. The above Littlewood-Paley decomposition is almost orthogonal in L^. 

Proposition 2.1. For any f G 5'(M^) and g G 5'(IR'^), the following properties hold: 

ApAJ = if \p-q\>2, 

Ag{Sp.ifApg)=0 if \p-q\>5. 

Having defined the linear operators Aq{q > —1), we give the definition of Besov spaces and 
Bony's decomposition. 

Definition 2.1. Let 1 < p < oo and s G M. For 1 <r < oo, Besov spaces C S' are defined 
by 

f G B;^, ^ ||/||b|,„ =: ( ' < oo 

q>-l 

and Bp ,^ C S' are defined by 



f G B;^oo ^ =: sup 2'1'\\AJ\\l, < oo. 

q>-l 



5 



Definition 2.2. Let f,g be two temperate distributions. The product f ■ g has the Bony's 
decomposition: 

f-g = Tfg + Tgf + Rif,g), 
where Tfg is paraproduct of g by f, 

Tfg= J2 ApfAgg = J2Sg-ifAgV 

p<q-2 q 

and the remainder R{f, g) is denoted by 

R{f,g) = ^^qf^q9 with := Ag_i + Ag + Ag+i. 

As regards the remainder of para-product, we have the following results. 
Proposition 2.2. Let (si,S2) € M? and 1 < p,pi,P2,r,ri,r2 < oo. Assume that 

111,111 

-< \ <1, -< 1 , and si + S2>0. 

p Pi P2 r ri r2 



Then the remainder R maps B5,^ x Bt,^ in Bp^r " ""^ and there exists a constant C 



Sl+S2+(i{i — 1 ^) 

777 h 

such that 

(j\si+S2\+l 



Some conclusions will be used in subsequent analysis. The first one is the classical Bernstein's 
inequality. 

Lemma 2.1. Let A; G N and < i?i < i?2- There exists a constant C, depending only on Ri,R2 
and d, such that for all 1 < a < b < oo and f G L", 

Supp JT/ c B(0,i?iA) ^ sup Wd^fhb < C'=+U'=+<i-i)||/||ia; 

|Q;|=ifc 

Supp Tf c C{0,RiX,R2X) => C-''-^X''\\f\\La < sup Wd'^fU-^ < C*^+^A'=||/||lc. 

|a|=A; 

Here T f represents the Fourier transform on f. 

As a direct corollary of the above inequality, we have 
Remark 2.1. For all multi-index a, it holds that 

I|9"/||b^,.<C||/||^.,h. 
The second one is the embedding properties in Besov spaces. 
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Lemma 2.2. Let s G M and I < p,r < oo, then 

Bt, ^ ^ Bf, whenever s < s or s = s and r < r: 

Bp J. ^ Bp^ ^ ^' whenever p > p; 

B^/f{l<p<oo)^Co, B^^.^CDL^, 
where Co is the space of continuous bounded functions which decay at infinity. 

The third one is the traditional Moser-type inequality. 

Proposition 2.3. Let s > and 1 < p,r < oo. Then Bp ^. L°° is an algebra. Furthermore, it 
holds that 

< C(||/||Loo||5r||Bs_^ + ||5||L-||/||iJ|,J. 

On the other hand, we present the definition of Chemin-Lerner's spaces first introduced by 
J.-Y. Chemin and N. Lerner [6j, which is the refinement of the spaces Lj,[Bp j.). 

Definition 2.3. For T > 0, s € M, 1 < r, p < oo, set (with the usual convention if r = oo) 

q>~l 

Then we define the space Lf^{Bp^^) as the completion of S over (0, T) X M'^ by the above norm. 
Furthermore, we define 

Ct{BI,) ■.= L^{Bl,)nC{[Q,T],Bl,), 

where the index T will be omitted when T = +oo. Let us emphasize that 
Remark 2.2. According to Minkowski's inequality, it holds that 

Then, we state the property of continuity for product in Chemin-Lerner's spaces Lj,{Bp^). 
Proposition 2.4. The following estimate holds: 

whenever s > 0,1 < p < oo, I < p, pi, P2, Ps, Pi ^ oo and 

11111 
- = — + — = — + —. 

P Pi P2 P3 Pi 

As a direct corollary, one has 

ll/5llzp(B=_,.) < '^ll/llz?,Mi?^.jll5llz?,2(B|,,) 
whenever s > N/p, - = — + —. 



In addition, the estimates of commutators in L^(i?p ^) spaces are also frequently used in the 
subsequent analysis. The indices s,p behave just as in the stationary case [U [12] whereas the 
time exponent p behaves according to Holder inequality. 

Lemma 2.3. Let 1 < p < oo and 1 < p < oo, then the following inequalities are true: 

' ^''i\\fh'^HB;J\9\\L'^HB;^,r s = l + N/p, 

< < CCq\\f\\r^n^^,^Jg\\^P2^j^s+ly s = N/p, 

^ C'Cg||/||jPi(^.+i)|bllzP2(5._^), s = N/p, 

where the commutator [•, •] is defined by [f,g] = fg — gf, the operator A = div or V, C is a 
generic constant, and Cq denotes a sequence such that ||(cg)||;i < 1, ^ = ^ + ^• 

In the symmetrization, we shall face with some composition functions. To estimate them, 
the following continuity result for compositions is necessary. 

Proposition 2.5. Let s > 0, 1 < p,r, p < oo, F e M) with F{0) = 0, T £ (0, oo] 

and V G L^{B^^^) n L^(L°^). Then 



\s] + l\ 



Finally, we give the estimate of heat equation to end up this section. 

~ s— 2+ — 

Proposition 2.6. Let s e R and 1 < a,p,r < oo. Let T > 0,-0,0 G B^,r ^"-^ / ^ L^{Bp^r 
Then the problem of heat equation 

dtu - p.Au = /, ^14=0 = uq 

has a unique solution u € L^{Bp^r'^) H L^{Bp j.) and there exists a constant C depending only 
on N and such that for all ai £ [a, +oo], we have 



i_ 

In addition, if r is finite then u belongs to C{[0,T];Bp j.). 



3 Well-posedness for k 

In this section, using the frequency-localization methods, we give the proof of main result. 
The proof of Theorem \l.l\ The coefficients Tp,Tyj,\ are assumed to be one. For classical solu- 
tions, (jl.ip can be changed into the following system in (n, u, T, 

dtu + div(nu) = 0, 

ndtu + n{u ■ V)u + V (uT) = nV ^ — nu, 

ndtT + nu • Vr + (7 - l)nrdivu = (7 - l)nAT + 2^n|u|2 - n{J - Tl), ^ ^ 

A <P = n — n. 
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In order to obtain the effective frequency-localization estimate on (n, u), we introduce a function 
change 



/ p \ 


/ In n — In n \ 


u 


u 


e 


T-n 


V E / 





Then the new variable {p,u,9,E) satisfies 

dtp + u ■ V/5 + divu = 0, 
dtu + TlVp + (u • V)u + V6 + eVp = E - u, 
dtO - ^^^Ae + u-Ve = hi{p)M - (7 - 1){Tl + ^)divu + 
t dtE = -VA-Miv(/i2(/5)u + nu), 



^lu|2 
2 l"l 



(3.2) 



where hi{p),h2{p) defined by 



hi{p) 



(7-1)k 



n 



1 — exp(— p) ) and /i2(/o) = n(exp(/9) — 1) 



are two smooth functions on the interval (—00, 00). The non-local term VA ^V-/ is the product 
of Riesz transforms on /. Here and below, we set h = ^"^"^^^^ for simplicity. 
The initial data (11.21) become 



(/f), u,6l,E)(x,0) = (In no - lnn,uo,7o - Tl,^A ^{riQ - n)), x G 



Remark 3.1. The variable transform is from the open set {(n, u,T, E) € (0, -|-oo) x IR" x 

X X M X 



(3.3) 

X 

R^} to the whole space {{p, u, 6*, E) € M x x M x R^}. It is easy to show that for classical 
solutions (n, u, T, E) away from vacuum, ()l.ip - (ll.2p is equivalent to (I3.2p - ()3.3p . 

The proof of the local well-posedness stems from a standard iterative process. First of all, 
we consider the linear coupled system of hyperbolic-parabolic form 

dtp + V • Vp + divu = 0, 

9tu + rLVp+(v v)u = /, 

dtO - RM = g, 
dtE = - VA-^div/i, 



(3.4) 



subject to the initial data 



(p,u,e,E)| 



i=0 



where y^, f,h 



pN 



and g : 



(Poj'^iOi^cEo) 

s>N .ire 



For the system ()3.4p - ()3.5p . we have the following conclusion. 
Proposition 3.1. Letp € [1, +00], r € [1, +00), si > and S2 S 



(3.5) 



Suppose that (/Jo,uo,Eo) € 



Sp^j.) and 



Vv € < 



' C{[0,T],B^\-^) if si > l + N/2,or si = l + N/2 and r = 1; 
C{[0,T],B2% ) for some e > if si = 1 + N/2 and r > 1; 

N_ 

. C([0,r],i?2%nL-) z/ 0<si<l + iV/2; 
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for any given T > 0. Then the system ^3.4\l-(375\) has a unique solution (p,u,6,'E) satisfying 



(p,u,E) G CriB'^^r) and 9 G CriB'p}.). 

Proof. Note that the L^- boundedness of Riesz transform, Proposition 13.11 is the direct conse- 
quence of Proposition 12.61 and Theorem 4.15 in the recent book [5]. □ 

In what follows, the proof of Theorem 1 1.1 1 is divided into several steps, since it is a bit longer. 
Stepl: approximate solutions 

We use a standard iterative process to build a solution. Starting from (p'', u'', 0'', E'') := 
(0,0,0,0). Then we define by induction a solution sequence {(p™, u"^, 0"^, E'")}mGN by solving 
the following linear equations 

f (9t/o™+^ + • Vp'"+^ + divu™+i = 0, 
(9tu™+i + TlV/9™+i + (u*" • V)u"^+i = - 6i™V/9™ + E™ - u™, 

= -u™ • V6i"^ + /ii(p™)A6l" - (7 - 1)(Tl + 6i'")divu™ + 2^|u™|2 - 6*™, 
I ^tE^+i = -VA-idiv{/i2(p"')u'^ + nu™}, 

with the initial data 

(p'™+\u"+\0'"+\E"^+!)(x,O) = (5™+ipo, 5^+iuo, 5^+1^0, WEo), x G M^. (3.7) 

Since all the data belong to B^^, Proposition 13.11 enable us to show by induction that the 
above Cauchy problem has a global solution which belongs to C{B^j.). 

Step2: uniform bounds 

Set 

:= Ct{BI,) X {CT{Bl,)y X Ct{B^+') X {cABl,))^ 

for r > 0. We hope to find a time T such that the approximate solution {{p"^, u™, 0™, E™)}mGN 
is uniformly bounded in E^. 

First, by applying the operator Ag(q > —1) to the first two equations of ()3.6p . we infer that 
for (Agp™+i, Agu"+!) 

' dtAqp""-^^ + (U"^ • V)Aqp"'+^ + A^divU^+l = [U™, Ag] • V/9'"+\ 

Aqu™+i + TlA^Vp^+i + (u™ • V)Agu'"+i (3-8) 



AqV6''" + [u™, Ag] • Vu"+i - Vp^A^ei"^ + [V/)'", Aq]e"' + AgE™ - A^u"^, 



where the commutator [•, •] is defined by [f,g] = fg — gf ■ 

Then multiplying the first equation of Eqs. (13. Sp by TLAgp™"*"^, the second one by Aqu'"^^, 
and adding the resulting equations together, after integrating it over R^, we have 



1 d 
2'dt 



^ " (t'.IIA I II a ii™+1||2 
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= i y divu"^(rL|Agp"^+il2 + |A,u"^+i|2) + I TlK, Ag] • Vp'"+^A,/5"^+i 

- J AgVe™ • Agu'"+i + J [u™, A,] • Vu'"+iA,u™+i - J Vp" • AgU™+i Ag^"^ 
+ J [Vp"", Aqje"" ■ Agu"+i + J AgE"^ • Agu"^+^ - y Agu*" • Agu'"+i 

< i||Vu"b^(||A,p™+i||i, + |lA,u"+ii,) + rL||[u™,A,].Vp™+ikH|A,p"+^^^ 

+ ||A,V0-|U.||A,u-+i||^2 + ||K,A,] • Vu-+i||i2||A,u-+i||i2 

+ ||Vp™||Loo||A,0-||^2||A,u"+i|Ii2 + ||[V/,-, A,]0™||i.||A,u-+i||i2 

+ ||A,E'"||^.||A,u'"+i||^2 + ||A,u-||i.||A,u"+i||^., (3.9) 

where we have used Cauchy-Schwartz's inequahty. 

1 

Dividing 1^ by (^Tl|| Agp^+i + ||AqU™+i||2 2 + ' (e > is a small quantity), we get 
d 



< C||Vu™||ioc(||A,/5-+i||i2 + \\AgU^+'\\L2) + C||[u-, A,] • Vp™+i| 



L2 



+C||A,V0™||i2 +C||[u", A,] • Vu-+i||i2 +C7||V/5"^||Loo||A,0™||i2 



V 1/ 11^^ -r ^ II , 



+C\\[Vp^,A,]e^h2 +C\\A,E^h2 +C\\A,u"'h2, (3.10) 

where C > here and below denotes a uniform constant independent of m. Integrating p.lOp 
with respect to the variable t £ [0,T], then taking e — > 0, we arrive at 

l|A,/."^+^(t)k2 + ||A,u-+i(t)||^. 

< C7(||A,p-+ii. + \\A,u^+'h2) + cJ^ ||Vu"(T)b^(||A,/,™+i(r)||^. 

+ \\A,n^+\T)h2)dT + C 1^ (l|[u",A,] ■ Vp'^+ii^ + ||[u-,A,] • Vu™+i||i.)dr 
+Cy* (||A,V0™||^2 + llVp-lli^llAge^lli^ + \\[Vp^,A,]e^h2 

+ ||A,E™||i2 + ||A,u'"||i2)dr. (3.11) 

Multiply the factor 29'^(cr = 1 + iV/2) on both sides of ([311]) to obtain 
2^'^||A,p™+i(t)||^2+2^'^||A,u'"+i(t)||i2 

< C2^'^(|lA,p™+i|l^. + IIA^ur^L^ + Cl^ ||u-(r)b.^2'"^(||A,;.-+i(r)||^. 

+C I' C,{T)\\p"'\\Bl,\\9"'\\Bl,dT + C j\'"'(\\A,Ve'^\\L2 

+ ||Vp"||Loo||Age'"||i2 + \\Aq-E^\\L2 + ||Agu"^|U2)dT (3.12) 
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where we used Remark 12. II and Lemmas 12.2112.31 and {cq} denotes some sequence which satisfies 
II (cq) 11(1 < 1 although each {cq} is possibly different in (I3.12|) . 
Summing up (|3.12p on g > — 1 implies 



ll(p"+\u'"+i; 



< C(||(po"^+i,uri)||^.^ ||u™(t)||B,^J|(/,™+\u™+i)(t)||j^^(^.^). 

+c£ ((1 + ||p"b,^,J||r'||^.+i + ||(u-,E-)|b.^)dt. (3.13) 



Then it follows from Gronwall's inequality that 



< Ce^^-(^){||(,o,uo)||Bj, 

a + IIP™WbfJII^"^WIIi,,^|^ + ll(u'",E-)(t)b.^)dt}, (3.14) 







with Z"(r) :=/; ||u-(t)||B./t. 

On the other hand, by the last equation of (|3.6p . with the aid of Proposition 12.51 we can 
obtain 

||Tpm+l 11 

< c(||E™+1||b.^+^ ||VA-idiv(/i2(p'")u- + nu-)||Bj/t) 

< c(^\\Bo\\bi,+ j\l + \\p"'\\Bi,)\\un\Bi,dty (3.15) 

where we have used the L^-boundedness of nonlocal (but zero order) operator VA~^div. 

Taking ai = a = oo, s = a + 1, p = 2 and r = 1 in Proposition 12.61 applying the 
resulting inequality to the third equation of (13. 6p . we have 



ll^^^^llz-CB-ji) ^ '^(IIC^^IIbJ+i + i^ + T)R ^\\K\\l^(B^-i)), (3-16) 

where := -u"^ • VO"" + hi{p'^)M'^ - (t - 1)(Tl + 6'"^)divu™ + - 0™. From 

Propositions 12. 4112. 51 and Lemma 12.21 we are led to 



gm+l I 



< C{\m^... + (1 + T)K-i((l + ||(p"^,u-)||£^(^.^p||0-||£^(^...^ 

+(l + ||u'"||j^(^.^p||u'"||j^(^.^))}. (3.17) 

Note that although the above constant C maybe depend on N , it is nothing to do with m, so 
we obtain the following uniform estimates. 
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Lemma 3.1. There exists a time Ti > (independent of m) such that 



(3.18) 



for all m € N U {0}, provided that k > is sufficiently large, where the constant Ci > 
independent of m and A := ||(/9oiUo,Eo)||_b^, + ||6'o|1 r<t+i . 

2,1 -^2,1 



Proof Indeed, the claim follows from the standard induction. First, we see that ()3.18p holds for 
m = 0. Suppose that (IS.lSp holds for any m > 0, we expect to prove it is also true for m + 1. 
Together with the assumption, by (|3.14p - (|3.15p and (|3.17p . we get 



1(A 



m+l ,,ni+l\|i ^ ^ 



A + Tmax{2CiA,2{CiAf} 



CCiAT 



(3.19) 



IE' 



m+l I 



'iT(Sf.l) 



< c 



A + Tmax{2CiA, 2{CiAf} 



(3.20) 



gm+l I 



< c 



A + {1 + T)k~^ max{4C7i A, 4(Ci A)^} 
A + rmax{4C7iA, 4.{CiAf} 



(3.21) 



where we suffice to take h satisfying k > ^^^^ (T to be determined) in the last step of the 
inequality (|3.2ip . Combining with (|3.19p - (|3.2ip . we have 



~ _l_ ll~ 



< 3C 



A + Tmax{4C7iA, 4:{CiAf} 



(3.22) 



Furthermore, if we choose T\ satisfying 



< Ti < min{ ^"(g^~j^\ ro}(Ci > 1 + 6C), 



where Tq is the root of algebra equation 

then ||(p-+\u-+\E'"+i)||j;^(^^.^) + 
proof of the assertion. 



6Cmax{4,4Ci>l}t ' 



□m+l 11 



u < C\A is followed, which concludes the 



□ 



That is, we find a time Ti > (independent of m) such that the sequence {(p™', u™', d'^ 
E'")}meN is uniformly bounded in Ej^^. 

Step3: convergence 

Next, it will be shown that {(p'", u™-, 6''", E'")}„jgN is a Cauchy sequence in E'^^. 
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(3.23) 



Define 

j^m+l _ pTn+p+l _ pm+l — _ 

_ Qm+p+1 _ Qm+1 — £jm+P+l _ jim+l 

for any {m,p) G N^. 

Take the difference between the equation (13. 6p for the (?tt- + |) + l)-th step and the {m + l)-th 
step to give 

' dt5p"'+^ + u'"+P • V(5/o™+i + (5u™ • + div(5u™+i = 0, 

+/ii(p"^)A50™ - (7 - l)(7i + 0™)div5u™ 
-(56i'"divu'™+P + l^(u™+P + u'")5u"^ - 56''", 

^ 9t5E"^+^ = -VA-Miv{[/i2(p'^+P) - /i2(p")]u™+P + /i2(p"')(5u'" + nJu™}, 

subject to the initial data 

(5p"+\5u'"+\5e"+\5E"+i)(x,0) = [5^+p+i-5™+i](po,uo,0o,Eo), x G M^. (3.24) 

Applying the operator Ag(g > —1) to the first two equations of p.23p gives 

(u'"+P . V)Ag(5/9™+i + Agdiv5u™+i 
V(5/9"^+i - Ag(5u™ • Vp'"+^), 

TlA,V(5p'"+i + (u™+P • V)A,(5u'"+^ (3-25) 
V^u'^+i - A5(5u"^ • Vu"^+i) - A^V^^"^ 
Aq{WVp"'+P) - V5p"A,6i™ + [V(5/9™, Ag]6i™ + A^JE™ - A,(5u"^, 

By multiplying the first equation of Eqs. p.25p by TLAq5p™+^, the second one by Aq5u'"+^, 
and adding the resulting equations together, after integrating it over M^, we have 



dt^qSp^- 


fl 


= [u"^+P, 


A 


StAgJu" 


+1 


= fu™+P, 


A 



■(rL||A,5p™+ii|i, + ||A,<5u"+i||i, 



1 d 
2Jt 

< \\Vu"'+P\\Lo.[TL\\AgSp^-^'\\l, + \\A,5u'^+%, 

+ (rL||[u™+f, A,] • V5/;™+i||i2 + \\A,{Su^ ■ Vp^+')h2)\\A,6p^+'\y 
A,] • V5u™+i||i2 + ||A,(<5u'" • Vu"+i)||i2 + \\AgV6e'^\\L2 
+ \\Agi6e'^Vp"'+P)\\L^ + ||5p'"||Loo(||A,Vr^||i2 +25||A5e-||^2) 
+ \\[V6p"',A^]e"^h2 + ||A,5E'"||i2 + ||A,<5u"||i2}||A,5u"+i||i2, (3.26) 
where we have bounded the integration 
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< ||<5p"||ioo(||A,V0'"||i2+25||A,e"|1^2)||A,,5u-+i||i2. 
Similar to the estimate of (|3.10|) . we can obtain 

|(rLi|A,5p-+i|li. + l|A,<^u-+^!li.+e)^ 

< !|Vu"+PiUoo(rL||A,5p'"+i||i2 + ||A,Ju'"+i||^2) + (Tl II [u™+P, A,] • V5p^+^\\l2 

+ ||A,(5u- • Vp"^+i)||i2) + {||[u'"+^ A,] • V5u-+i||i2 + ||A,(5u- • Vu"+i)||i2 

+ \WV5e^\\L2 + \W{5ff^Vp^+n\\L2 + ||5p"||Loo(||A,V0™||i2 +2«||A,e"||i2) 

+ ||[V5/9-,A,]e-||i2 + ||A,5E-||i2 + ||A,<5u"||i2}, (3.27) 

where e > is a small quantity. 

Integrating (|3.27p with respect to the variable t G [0, Ti], then taking e — >■ 0, we arrive at 

< c(||A,5p™+i^. + ||A,<^u-+i||^.)+C^Vu'"+niLoc(||A,5p-+i||^. 

+ ||A,5u-+i||i2) +C^* (||[u-+P,A,] ■ V<5p-+i||,.2 + ||A,(<5u- • Vp-+i)||^2) 

+C y* (||[u'^+P,A,] • V^u'^+illia + ||Ag(5u™- Vu™+i)||i2 + ||A,V50™||i2 

+ ||A,((50™Vp-+P)||^2 + ||5p™||ioo(||A,Ve"'||i2 +2^||A,e-||i2) 

+ ||[V5p™, A,]^™||^2 + ||A,5E"^||i2 + ||A,<5u™||i2). (3.28) 

By multiplying the factor 2'^^°'^^^ on both sides of the resulting inequality (j3.28p . we obtain 
2^(^-1) (||A,5p-+i(t)||i2 + \\A,6u"^+\t)h.) 

< C2^('^-i)(||A,<5priL2 + IIA^c^u^+Il^) 

+cl^ ||Vu-+p||lo.2^('^-i)(||A,5p™+1||^2 + \\A,Sn^+'h2) 

+C f (c,||u-+^b.J|5p™+i^^_, +c,Pu-||^.-.||p-+i||,^. J 

+cj|'^^'"b,^,, +c,||<^^™b,^J|p^'"+niB,^,, +c,||<^P-||l^(||V0-||^.^ + ii^'-b,^,,) 
+cJ|V^'"||^.-iPp™||^.-i +2^('^-i)||A,<5E™||^2 +2''('^-i)||A,5u™||i2), (3.29) 
where {cq} denotes some sequence which satisfies ||(cg)||ii < 1. 
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Summing up (|3.29p on g > —1, it is not difficult to get 



||(5p'"+^5u-+l)||r.c.o.-l^ 



Jo 



< C72-™|l(/,o,uo)b,^^^ 



' 

X (l + II , p^+P, u-+^ , e™) II ^it, (3.30) 



where we liave used Lemma 12.11 and Remark 12.11 

Witli tiie aid of Gronwall's inequality, we immediately deduce that 



< Ce^^-''m){2-||(^o,uo)bj^^ 



rp 

+ I 'e-^^'"''''(*)(||(5p'",<^u-,,5E-)||^.-i + ||50-b,^^ 



x(l + ||(/,-+\/,-+P,u-+\nb?,,)rfi} 
< Ce^^^{2-- + T,[\\{6p-,6u"^JB"^)\\~^^^^^^^^^^^ (3.31) 

where we have noticed Remark 12.21 and the fact that the sequence {(p™, u*", 0*", E™)}mgN is 
uniformly bounded in E^^ . 

From the last equation of (|3.23p . we get directly 



< mrx.-^ + / ' (ii5p"iuj^^i||u-+^iib,^,, + ii5u™ii^.-.(i + \\pn\Bi,] 

< C72-™ + CTi||(5p",5u'")||;^^(^.^i). (3.32) 
Using Proposition 12.61 (taking ai = a = oo, s = a, p = 2 and r = 1), we have 

II'5^"^^IIl-(bj,) 
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< c[\\6e^^'\\B^^^^ + (1 + ri)/i-i||Fr||jo.(5.^2)), (3.33) 

where 

(/>'") A J0™ - (7 - 1){Tl + e"')div5u" - Je^divu'^+P 

+l^{u'^+p + u"')5u" - se"". 

In bounding F^, each product term can be estimated effectively with the help of the standard 
Moser-type inequality (Proposition 12. 3( ). except for the term hi{p"'')A69^\ Here, we develop a 
Moser-type inequality of general form to estimate hi{p"^)A69"^, which will be shown in the 
Appendix, see Proposition 14.11 According to it, we can reach 

||/ii(/9'")A5e"^||^.-2 < C(||/ii(p")||loo||A(50'™||^.-2 + ||A50™||i2||/ii(/9"^)|l^^-2). (3.34) 
The second term in (I3.34p can be further estimated as 

||A50-||i2||/ll(p™)||^.-2 

00,1 

< C7||A(^0™||^.-2||/ii(p™)||^A.-i {a = l + N/2, N > 2) 

< c\\5e"'\\Bi,\\hi{p'^)\\Bi, 

< C\\6e'''\\Bi^\\p"'\\Bi,, ' (3.35) 

where we used the embedding properties -62!^ ^ ^^'^ ^21^ ^ -^00 i- ensure ^ 
B^{^, - 1 < 1 + N/2 i.e. A^ < 4 is required in the last second step of 

Combining (j3.33p - (|3.35p and recalling on the choice of «;(k > ^^yp'), we conclude that 

II r/im+l 1 1 

< c{2~- + T,{\\{5p-Mn\\l^^^B^^-^^ + ¥^^^^ (3.36) 
Therefore, together with (|3.31|) - (|3.32|) and (|3.36p . we end up with 

|[(5p-+\5u'"+\5E™+i)||j^^(^.-.) + l|5^™+iz?.(B,^,) 
< CT\2-- + T,(\\{5p-M'^,5^n\\l^^^B^^-^^^ (3.37) 
where Ct^ ■= Ce^^^. Arguing by induction, one can easily deduce that 

||(<5p-+\<5u™+\5E™+i)||^^^(^.-,) + \\S9-^%^^^Bi,) 

+Ct,^2-('"-'=)^^1^. (3.38) 



fc=0 
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As \\{5pP ,5u'P,5'EP)\\joc ( p,cr-i-. + \\S9P\\jao (f>a- ) Can bc bounded independent of p, We fuither takc 



Ti SO small that 

< C2-". a.d ^^-'^f"'' < 4-'. 

(m + 1)! ~ A;! ~ 

Thus we conclude that there exists some constant C2 > (independent of m) such that 

||(Jp-+\<5u™+\<5E™+i)||£^^(^._.) + < (3.39) 

which implies {(/>™, u™, 0™, E™)} mgp^ is a Cauchy sequence in E!^ ^. Therefore, there exists 
some function (p, u,0,E) in E!^~^ such that 

{(/o"',u™,0™,E™)} ^ (p,u,0,E) strongly in 

Step4: the solution [p, u, 6, E) 

In this step we show that (p, u, 6',E) G E^_^ is a solution of the system (j3.2p - p.3p . Fatou's 
property ensures that {p,u,9,B) G L^^{B^^^) x (Z~(5f J)^ x L5?(50^) x (L~(5f J)^, since 
{(p"',u™,6'"^,E™)}^eN is also uniformly bounded in the spaces L??(B^,i) x (L??(B^,i))^ x 
Z~(50i)x(Z~(52-,i))^. 

On the other hand, {{p"^, u™, E™)}„gN converges to (p, u, E) in C([0, Ti]; B^;^ ^) and {O'^jmeN 
converges to 6 in C([0, Ti]; i?2 i)- These properties of strong convergence enable us to pass to 
the limits in the system (j3.2p - (j3.3p and conclude that (p, u, 0,E) to the system (|3.2p - (j3.3p . 
Now, what remains is to check {p,u,e,E) also belongs to C([0, Ti]; J x (C([0, Ti]; ^^J)^ x 
C([0, Ti]; i?2 i"^) X (C([0, Ti]; i?2 i))'^- Indeed, for instance, we easily achieve that the map 
t I-)- II Aqp(t) |1j;^2 is continuous on [0,Ti], since p € C([0, Ti]; i?2l^)- Then we have Agp(t) G 
C([0,Ti];Sf J for all q > -1. Note that p G Z5?(Sf i), the series Eg>-i 2'^'^ 11^9/0(0 IIl2 con- 
verges uniformly on [0,Ti], which yields p G C([0, Ti]; i?2 i)- The same arguments are valid for 
the other variables (u,0,E). Hence, we finish the existence part of solutions. 

Step5: uniqueness 

Let p = pi — P2, u = ui — U2, = 9i — 62, E = El — E2 where (/Oi, ui, 6*1, Ei)^ and 
(p2; U2; ^2, E2)^ are two solutions to the system ()3.2p - ()3.3p subject to the same initial data, 
respectively. Then the error solution (p, u, ^,E)^ satisfies 

dtp + divu = — ui • Vp — u • Vp2, 

atu + TlVp = -V^- ui • Vu - uVu2 - 0iVp- eVp2 + E - u, 

dt9-KAe = -ui • Ve- uVe2 + [hiipi) - hi{p2)]Aei + hi{p2)M (3.40) 
-(7 - 1)(Tl + ^i)div5 - (7 - l)Mivu2 + ^5(ui + us) - 6, 

_ 5iE = -VA-iy • [(/i2(pi) - /i2(p2))ui + (/i2(/>2) + n)S]. 

As previously, following from the proof of Cauchy sequence, we obtain the inequalities: 

ll(?'")IIZ5?(B-7i) ^ ^ (ll(P>U,E)||^.-i) (^1 + ||(pi,/32,Ui,U2,6li)||B-^)dt 

+\mLi,^iBi,) (3.41) 
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and 

/•Ti 

IE 

'0 



Zoo(B-i) ^^l' IKp, u) 11^.-1(1 + II (p2,ui)||Bj^,)di. (3.42) 
According to Proposition 12.61 (taking ai = Q = l, s = a — 2, p = 2 and r = 1), we have 

\\^h},^iBl,) ^ ll(P2,Ul,U2)||joo(5.-i))||^||£i^^(5.^) 

+ ^^^|(p,5)||^.-i(1 + ||(ui,U2)||b^.^^ + 11(^1, 02)lbj,^^ (3.43) 



We further choose Ti satisfying 



2C(1 + ||(p2,Ui,U2)||j«.(^.-l))' 

and obtain 

ll^llzi (BM ^ C r'||(p,5)|| . 1(1 + ||(ui,U2)|b,^ + 11(01, ^2)bjjdt. (3.44) 

Ti^ 2,1^ 2,1 

Note that L^^{B^^^) = L\^{B^^^), we insert ([3:11) into ([OT]) to get after combining ([3:12]) 



< c 



I ' ll(p,u,E)||£^^(^.-i^(l + ||(pi,P2,ui,U2)||bj_^ + \\{ei,e2)\\B.+i)dt. (3.45) 



Gronwah's inequality gives (/), u, E) = immediately. Substituting it into ([3.44p . = is also 
followed. 

Finally, by Remark 13. 11 we can arrive at Theorem [LT] satisfying the inequality ([1.4p . Hence 
the proof of Theorem 11.11 is complete. 



4 Appendix 



In the last section, we give the crucial Moser-type inequality in the non-homogeneous Besov 
spaces and Chemin-Lerner's spaces. For the homogeneous version, which has been remarked by 
Zhou in 1361. 



Proposition 4.1. Let s > and 1 < p, r, pi,p2,P3,P4 < 00. Assume that f G L^^ PI Bp^^^ and 
g e LP^n B^^ ^ with 

1 _ 1 1 _ 1 1 
P Pi P2 P3 Pi' 

Then it holds that 

\\fg\\Bi^ < C(||/||lpi II^IIb^,,. + II5||lp3 ||/lb^^,J. (4.1) 

In particular, whenever s > N/p, there holds 

||/9||b|,. <C||/||B^ Jb||B|,.. (4.2) 
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Proof. From Bony's decomposition, we have 

fg = Tfg + Tgf + R{f,g). 
It follows from Proposition 12. II that 

2'^^\A,Tfg\\Lr < 2^^ \\MSg,fA,,g)\\Lv 



|g-g'|<4 

< V 2('?-«')^||S,,/||z.«2^'«||A,,5||lp2 (- = - + - 
|,-7t<4 P^' 

< C J2 2(^-'?')l|/||iPi2'?'^||A,,5l|LP2 

\q-q'\<4 

< Ccgi\\f\\L.i\\g\\Bs^, (4.3) 



where Cgi := Y.\q-q'\<4—§^^^^^ satisfies ||cgi||p- < 1. Similarly, 



2'''\\A,Tgf\\LP < 2'?^ Yl \\^q{Sq'9\'mLp 

\q-q'\<i 

< Cc,2\\9\Mf\\Bs„., (4.4) 



where Cg2 := E|g-n'|<4 9|| fy"!^"^"^ satisfies ||cg2||£'- < 1- 
On the other hand, from Proposition 12.21 we arrive at 

2'^'\\AgR{f,g)\\Lr 

< Cc,3p(/,5)b|„ 

< CcMso^Jg\\Bs^Js>0) 



< Cc^sll/llLPilblbg,,, (4.5) 



P2'1 



where Cqs := — ||i?.(/ 9)\\b^ ^^^^^^^^ llcgsll^'- < 1- In the last step, we used the embedding 
property B^_^ i L^^ 



Pi,oo- 



Hence, (j4.ip follows from (|4.3p - (|4.5p . Moreover, the embedding properties in Lemma [2. 2 1 give 
2]) directly. □ 



It is not difficult to generalize Proposition 14.11 to the framework of Chemin-Lerner's spaces 
Lj,{Bp,^). The indices s,p,r behave just as the stationary case whereas the time exponent p 
behaves according to Holder's inequality, which is given by a proposition for clarity. 

Proposition 4.2. The following estimate holds: 

II/5||ZP(B._J < C(||/||iPi(iPi)|bllzP2(B.^ + llffllL?,3(LP3)ll/llz?4(B.^ ,,)) 

whenever s > 0, l<p,r<oo, 1 < Pi,P2,P3,P4 1^ oo, 1 < p, pi, P2, P3, P4 < with 

1 _ 1 1 _ 1 1 
P Pi P2 P3 Pi 
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and 

1 _ 1 1 _ 1 1 
P Pi P2 P3 Pa' 

As a direct corollary, one has 
whenever s > N/p. 
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